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Abstract. Due to the finite size effects, the localisation of the phase transition in finite systems and 
the determination of its order, become an extremely difficult task, even in the simplest known cases. In 
order to identify and locate the finite volume transition point To{V) of the QCD deconfinement phase 
transition to a Colorless QGP, we have developed a new approach using the finite size cumulant expansion 
of the order parameter and the Lmn-method. The first six cumulants (71,2,3,4,5,6 with the corresponding 
under-normalized ratios(skewness S, kurtosis k ,pentosis 77± and hexosis Hi, 2 , 3 ) and three unnormalized 
combinations of them {O = U = J\f = (t^k) are calculated and studied as functions 

of {T,V). A new approach, unifying in a clear and consistent way the definitions of cumulant ratios, is 
proposed. A numerical FSS analysis of the obtained results has allowed us to locate accurately the finite 
volume transition point. The extracted transition temperature value To{V) agrees with that expected 
(V) from the order parameter and the thermal susceptibility xt iT,V), according to the standard 
procedure of localization to within about 2%. In addition to this, a very good correlation factor is obtained 
proving the validity of our cumulants method. The agreement of our results with those obtained by means 
of other models is remarkable. 

PACS. 1 2.38.Mh,12.38.Aw,25.75.Nq,64.60.an 


1 Introduction 

1.1 Phase Transitions and Finite Size Scaling(FSS) 

During the evolution of our beautiful universe from the 
big-bang instant until now many phase transitions have 
occurred at different space-time scales.For this reason, the 
physics of phase transitions phenomena is considered in 
general to be a subject of great interest to physicists. It is 
easy to understand the importance of this subject because 
firstly, the list of systems exhibiting interesting phase tran¬ 
sitions continues to expand, including the Universe itself, 
and secondly the theoretical framework of equilibrium sta¬ 
tistical mechanics has found applications in very different 
areas of physics like string field theories, cosmology, ele¬ 
mentary particle physics, physics of the chaos, condensed 
matter ... etc. Phase transitions occur in nature in a great 
variety of systems and under a very wide range of condi¬ 
tions. 

Phase transitions are abrupt changes in the global behav¬ 
ior and in the qualitative properties of a system when 
certain parameters pass through particular values. At the 
transition point, the system exhibits, by definition, a sin¬ 
gular behavior. As one passes through the transition re¬ 
gion the system moves between analytically distinct parts 


of the phase diagram. Depending on which external pa¬ 
rameter of interest, there are various measurable quanti¬ 
ties which are based on the reaction of a system to its 
change. We call them Response Functions (RF). If the 
external parameter corresponds to the temperature, then 
the response function is called Thermal Response Func¬ 
tion (TRF). Technically, temperature driven phase tran¬ 
sitions are characterized by the appearance of singularities 
in some TRF, only in the thermodynamic limit where the 
volume V and the number of particles N go to infinity, 
while the density p = N/V remains constant. That is, at 
the transition point, some global behavior is not analytic 
in the infinite volume limit. This singularity is according 
to the standard classification [T] given by the (5—function 
for a first-order phase transition, while for a continuous 
phase transition (second-order), the singularity has the 
form of a power-law function. We shall frequently refer 
to the concepts of transition region and transition point 
in the case of a first order phase transition. By against, 
in the case of a second order phase transition, we rather 
use the concept of critical region and critical point. The 
singularity in a first order phase transition is entirely due 
to the phase coexistence phenomenon, for against the di¬ 
vergence in a second-order phase transition is intimately 
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caused by the divergence of the correlation length. Now, 
if the volume is finite at least in one dimension with a 
characteristic size L — ^ the singularity is smeared 

out into a peak with finite mathematical properties and 
Four Finite Size Effects(4FSE) can be observed [5]: 

(1) the rounding effect of the discontinuities, 

(2) the smearing effect of the singularities, 

(3) the shifting effect of the transition point, 

(4) and the widening effect of the transition region around 
the transition point. 

These 4FSE have an important consequence putting the 
first and the second order phase transitions on an equal 
footing. The behavior of any physical quantity at the first- 
order phase transition is qualitatively similar to that of the 
second-order phase transition. However, even in such a sit¬ 
uation, it is possible to obtain information on the critical 
behavior. Large but finite systems show a universal be¬ 
havior called “Finite-Size Scaling” (FSS), allowing to put 
all the physical systems undergoing a phase transition in 
a certain number of universality classes. The systems in 
a given universality class display the same critical behav¬ 
ior, meaning that certain dimensionless quantities have 
the same values for all these systems. Critical exponents 
are an example of these universal quantities. The knowl¬ 
edge of the finite-size dependence of the various TRF in 
the vicinity of the phase transition region provides a very 
important way to compute, using finite size scaling extrap¬ 
olation, the properties of systems in the thermodynamic 
limit. 


1.2 Finite Size Effects(FSE) in QCD Deconfinement 
Phase Transition 

It is well established that Quantum Chromo-Dynamics 
(QCD) at finite temperature exhibits a typical behavior of 
a system with a phase transition. At sufficiently high tem¬ 
peratures and/or densities, quarks and gluons are no more 
confined into hadrons, and strongly interacting matter 
seems to undergo a phase transition from hadronic state to 
what has been called the Quark Gluon Plasma(QGP) or 
’’Partonic Plasma” (PP). This is a logical consequence of 
the parton level of the matter’s structure and of the strong 
interactions dynamics described by the QCD theory [3]. 
The occurrence of this phase transition is important from 
a conceptual point of view, as it implies the existence of a 
novel state of matter, believed present in the early universe 
up to times ~ 10“®s. Indeed, the only available experi¬ 
mental way to study this QCD phase transition is to try 
to create in a laboratory, using ultra-relativistic heavy-ion 
collisions(URHIC), conditions similar to those in the early 
moments of the universe, right after the Big Bang. Due to 
its similarity to the early universe, an URHIC is often 
referred to as ” little bang”. The analysis of the whole re¬ 
sults obtained in all experiments at SPS, RHIC and LHC 
revealed that indeed a new state of matter is formed, con¬ 
sisting of a strongly interacting partons [1] . The existence 
of this finite volume hot deconfined matter is strongly in¬ 
dicated because some important signatures are observed. 
One example is the jet quenching phenomenon. Accord¬ 


ing to QCD, high-momentum colored partons produced in 
the initial stage of a nucleus-nucleus collision will undergo 
multiple interactions inside the finite volume collision re¬ 
gion, generating a parton shower before hadronization. 
Due to thermal effects the cross section of the hadrons 
formation and the fragmentation process decrease [Slim 
[7] and to the color confinement property of QCD, only 
the color singlet part of the quark configurations would 
manifest themselves as physically observed particles. All 
hadrons created in the final stage are colorless. Therefore 
the whole partonic plasma fireball needs to be in a color 
singlet state called Colorless QGP (CQGP). For this rea¬ 
son, one can consider the QCD deconfinement phase tran¬ 
sition as a transition from local color confinement (d^lfm) 
to global color confinement(d;^lfm). Lattice QCD, a the¬ 
ory formulated on lattice of points in space and time, is 
an other important framework for investigation of non- 
perturbative phenomena such as confinement and decon¬ 
finement of partons, which are intractable by means of 
analytic quantum field theories. As is well known, the lat¬ 
tice’s space-time volume is finite. Whereby in both cases of 
experimental and lattice simulation models, we are deal¬ 
ing with finite systems and, therefore, they require the 
development of theoretical approaches that can rigorously 
define the phase transition in a finite volume taking into 
account the color singlet condition. Locating the finite vol¬ 
ume QCD transition point is a challenge in both theoret¬ 
ical and experimental physics. 

1.3 Motivation 

In the thermodynamic limit there is no problem to locate 
the transition point since it manifests itself as a singu¬ 
larity point. By cons, in finite volume this singularity is 
smoothed and is shifted away, consequently the location 
of the phase transition and the determination of its order 
become very difficult. The idea of a phase transition is al¬ 
ways related to the idea of locating the transition point. 
Two fundamental questions appear to be very important 
that we try to answer in the present work. Firstly, how to 
locate the transition point in finite systems? And secondly, 
how can we say for sure that a certain physical quantity 
has a particular behavior when approaching certain point, 
which may be conceived as the transition point? It is im¬ 
portant to have a precise knowledge of the region around 
the transition point since many quantities of physical in¬ 
terest are just defined in the vicinity of this point. It there¬ 
fore seems very important to find more sensible quantities 
to construct new definitions of the finite volume transi¬ 
tion point involving a minimum of corrections. Recently, 
many works have shown the importance of studying the 
high order cumulants of thermodynamic fluctuations. For 
this reason and even in the finite volume case higher-order 
cumulants and/or generalized ratios of them have been 
suggested as suitable quantities because they are highly 
related to the nature of the phase transition and serve 
as good indicators for a real location of the finite volume 
transition point. Mathematically speaking, the thermody¬ 
namical fluctuations of any quantity are quantified by cu- 
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mulants in statistics and are related to generalized ratios 
of them. Generally, they are defined as derivatives of the 
logarithm of the partition function with respect to the 
appropriate chemical potentials. The cumulant expansion 
method is then considered by many physicists to be very 
sensitive to the behavior of the system in the transition 
region and then is viewed as a promising powerful method 
to analyse the deconfinement phase transition in finite sys¬ 
tem [Miisn]- Therefore finding new observables to permit 
us an accurate localization of the transition point in QCD 
phase diagram is more than necessary. From our hadronic 
probability density function (hpdf) which is related to the 
total partition function and which contains the whole in¬ 
formation about the phase transition as pointed firstly by 
Gibbs [5] , it seems logical to believe that this information 
survives when the volume of the system becomes finite. 
Our basic postulate is that it should be possible to lo¬ 
cate the finite volume transition point by defining it as a 
particular point in each term of the finite size cumulant 
expansion of the order parameter, suggesting a new ap¬ 
proach to solve the problem. We believe that the finite 
volume cumulant expansion, should show some charac¬ 
teristics as signals of the finite volume transition point. 
Indeed and in order to identify and locate the finite vol¬ 
ume transition point Tq{V) of the QGD deconfinement 
phase transition, we have developed a new approach us¬ 
ing the finite size cumulant expansion of the order param¬ 
eter with the Lmn-method [5] whose definition has been 
slightly modified. The two main outcomes of the present 
work are: 1) the finite size cumulant expansion of our hpdf 
gives better estimations, than the Binder cumulant m, 
for the transition point and even for very small systems. 2) 
the singularity of the phase transition in thermodynamic 
limit survives in a clear way even when the volume of the 
system becomes finite. 


2 Statistical Description of the System 
Containing the Hadronic Phase and the 
Colorless QGP 

2.1 Exact Colorless Partition Function 

In our previous work, a new method was developed which 
has allowed us to accurately calculate physical quantities 
which describe efficiently the deconfinement phase transi¬ 
tion within the Colorless-MIT bag model using a mixed 
phase system evolving in a finite total volume V [2]. The 
fraction of volume (defined by the parameter h) occu¬ 
pied by the HG phase is given by : Vhg = hF, and 
the remaining volume: Vqgp = (1 ~ h.)V contains then 
the GQGP phase. To study the effects of volume finite¬ 
ness on the thermal deconfinement phase transition within 
the QCD model chosen, we will examine in the following 
the behavior of some TRF of the system at a vanishing 
chemical potential (^ = 0), considering the two lightest 
quarks u and d {Nf = 2), and using the common value 
= 145MeF for the bag constant. In the case of a 
non-interacting phases, the total partition function of the 


system can be written as follows: 

^TOT(h, F, T, /x) = ZcQGp{^)ZHGO^)Zvac{i^), (1) 
where, 

Zvacih, F T) = exp(-(l - h)BV/T), (2) 

accounts for the confinement of quarks and gluons by the 
real vacuum pressure exerted on the perturbative vacuum 
(B) of the bag model. For the HG phase, the partition 
function is just calculated for a pionic gas and is simply 
given by, 

ZHG{h, F, T) = exp aHGhVT\ (3) 

The exact partition function for a GQGP contained in 
a volume Vqgp, at temperature T and quark chemical 
potential /i, is determined by: 


Zgqgp{T,Vqgp, fJ') 


xTr 



— TT —TT 

/X (^Nq - ^ -h iipTs -f iipYs^ , 


where M (ip, ip) is the weight function (Haar measure) given 
by: 


M{ip,ip) = 


sin I l{,ip + ^)] sin(^)sin ( 


/3 = y (with the units chosen as: ks = h = c = 1), and 
Hq is the free quark-gluon Hamiltonian, Nq denotes 

the (anti-) quark number operator, and I 3 and Fs are the 
color “isospin” and “hypercharge” operators respectively. 
Its final expression, in the massless limit, can be put in 
the form, 


Zgqgp{T, Vqgp,p) 


4 




dipdipM[ip, Ip) 


X ^0{v,4>,ip)VQGpT^ ^ 


( 6 ) 


with, 

= 5(O,O,^)-k0QG(<P,-0,^)- (7) 

The two functions are given in terms of (T, F, p) variables 
as follows: 


rin n a , 
b(0,0, -) = aQG + 


7r2T2' 


( 8 ) 


and 


GQG{y^,tp, ^) 


24 ^ 

q=r,b,g 




- 1 


24 ^ 
9=1 



(9) 
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The two factors B-hg and slqg which are related to the 
degeneracies of the particles in the system are given by, 

f age = f2(^^Q + 

^ B-HG '§0^'^ 

dg = 2Nf , dG = 2. and djr = 5 being the degener¬ 
acy factors of quarks, gluons and pions respectively, aq 
{q = r, b, g) are the angles determined by the eigenvalues 
of the color charge operators in eq. ©: 


in lb ip lb 2ib 



( 11 ) 


and ag {g = 1, ...,4) being: ai = ar — ag, a 2 = Og — 
c^b, 0:3 = at — ttr, 04 = 0. Thus, the partition function of 
the CQGP is then given by. 


Zgqgp (h) = Zqop (q) ZcG (q), (12) 


where 


^Gc(q) = ^ 

/ +7r /*+'7r ^ 

-TT J —7T 

is the colorless part and, 

Zqgp (h) = exp (1 - h)VT"5QG(0, 0, ^). (14) 

is the QGP part without the colorless condition. Finally 
the exact total partition function with the colorless con¬ 
dition is given by, 


Since our hpdf is directly related to the partition func¬ 
tion of the system, it is believed that the whole infor¬ 
mation concerning the deconfinement phase transition is 
self-contained in this hpdf. This hpdf should certainly have 
different behavior in both sides of the phase transition and 
then we should be able to locate the transition point just 
by analyzing some of its basic properties. Then we can 
perform the calculation of the mean value of any thermo¬ 
dynamic quantity Q{T, /i, V) characterizing the system in 
the state h by, 

1 

(Q(T, M, V)) = Jq (h, T, g, V) p (h) dh. (19) 
0 

In our previous work, as mentioned above, a new method 
was developed, which has allowed us to calculate eas¬ 
ily physical quantities describing well the deconfinement 
phase transition to a CQGP in a finite volume V [2]. The 
most important result consists in the fact that practically 
all thermal response functions calculated in this context 
can be simply expressed as a function of only a certain 
double integral coefficient Lmn- The principal idea of these 
Lmn has emerged in the beginning when we performed the 
calculation of the < h(r, V) > and then we consider that 
it will be very interesting if we chose the definition of L^n 
in a judicious way so that all thermodynamic quantities 
can, in one way or the other, be written as function of 
these Lmn’s: 


L 


m,n 




dipd'ibM (ifi, ib){g{ip, - 0 , 0 ))™ 


gq R(v,V-;T.v) 


( 20 ) 


Ztot (h) = Zq (h) ZcG (h) 


(Ih) where the function R (ip, ib] T, V) is given by. 


with, 

Zo (q) = ZHG{^)Zvac{^)ZQGp{^)- (16) 

This latter is only the total partition function of the sys¬ 
tem without the colorless condition, which can be rewrit¬ 
ten in its most familiar form obtained in earliest papers 

0: 


LnZo(r, P, h) = 


-b 


f 2 M 




rji2 


B 

rp4 


|(1 - h) - Oi/ch 


id 

2 Ti 


VT^ (17) 


2.2 Finite Size Hadronic Probability Density Function 
and Lmn-Method 


R (:^, i; T, V) = (g{p, i, 0) - ghg - VT\ (21) 

We can clearly see that these Lm,n (q) can be considered 
as a state function depending on (T ,P) and of course 
on state variable q, and they can be calculated numeri¬ 
cally at each temperature T and volume V. As we will see 
later the mean value of any physical quantity Q(T, /r, V) 
can therefore be calculated as a simple function of these 
Lm,n (q) evaluated in the hadronic phase: Lm,n (0) and 
in the CQGP phase: Lm,n (!)• An other important prop¬ 
erty of these Lmn coefficients coefficients relies on the fact 
that any derivative within the T variable and V variable 
giving rise to other Lm,n (q) coefficients, it is like making 
a connection between different (q) and mixing them 
in a simple recurrent relation m- 


The definition of the Hadronic Probability Density Func¬ 
tion in our model is given by, 


p{h) 


Z{h) 

1 

JZ(h)dh 

0 


(18) 


2.3 Reminder of Some Thermal Response Functions 
obtained previously 

The first quantity of interest for our study was the mean 
value of the hadronic volume fraction < h(T, P) >, which 
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can be considered as the order parameter for the phase 
transition investigated in this work. According to (HID, 
< h(T, V) > has been expressed as pi fTUlfTT] : 


< h{T,V) >= 


Lo2 (1) — Lq2 (0) — Lqi (0) 

Lqi (1) — Lqi (0) 


( 22 ) 


which shows the two limiting behaviors when approaching 
the thermodynamical limit : 


lim < h(T, V) >= 0, lim < h(r, V) >= 1. 

(T)^cx. (T)^ 0 

(23) 

The asymptotic behaviors of < h(T, V) >, can be related 
analytically to the Heaviside step function in the thermo¬ 
dynamical limit : 


lim (h(T,V)) = l-0(T-ro(oo)). (24) 

(y)—)-oo 


The second quantity of interest was the energy density 
e(T, V), whose mean value was also calculated in the same 
way, and was found to be related to < h(r, V) > by the 
expression, 


As well known, the set of moments fully characterizes a 
probability density function provided that they are all fi¬ 
nite. At the same time the set of cumulants that is another 
alternative and, for some problems is a more convenient 
description. Once the set of moments are known, the prob¬ 
ability distribution may be obtained via reverse Fourier 
transform, that is the function 17 (t) which is nothing that 
the mean value of the depending only on the t variable 
and called the characteristic function of the distribution 

fix): 


n{t) 


+ 00 



— oo 


f{x)d{x) = 1 + XI ^ ■ 


(27) 

So, once Q{t) is known, all moments are known. New coef¬ 
ficients C„, which were introduced by Thiele , can 

be defined from the Maclaurin development of the In [2{t) 


oo p 

In f2(t) = X ^ ■ ( 28 ) 

n—1 


rp 2 

<e{T,V) >=Y< 


dLnZ \ 
dT ) 


> 


(25) 


From our FSS analysis of the whole results, the 4FSE have 
been observed M- These same effects have also been 
noticed in the present work. We also wish to recall the 
definitions of the specific heat ct (T, V) = ^^^^^;^and 

the thermal susceptibility xt iT,V) = repre¬ 

senting the thermal derivatives of both < e(r, V) > and 
< h(T, V) >. These TRF are very sensitive to the phase 
transition. 


They are called the semi-invariants or cumulants of the 
distribution f{x). In another way when we define the cen¬ 
tral moments Mn, relatively to the mean value of x (ai = 
(x) )as. 


M, 


+ 00 

= / (a:-ai) 


' r/ (—l)^n! fc 

f{x)dx = - 107 (^ 1 ) 


k^O 


k\ (n — k)\ 


(29) 

Using (EID, dMD and (1291) one can easily express the cumu¬ 
lants C„ and the central moments Mn via the moments 


3 Finite Size Cumulant Expansion: 

Theoretical Calculations 

3.1 Definitions of the Moments, Central Moments,and 
Cumulants 

Let us briefly recall the standard cumulant expansion and 
review some of its main properties.In probability theory 
and statistics, the cumulants Cn of a probability distribu¬ 
tion are a set of quantities that provide an alternative to 
the moments of the distribution. The moments determine 
the cumulants in the sense that, any two probability dis¬ 
tributions whose moments are identical, have identical cu¬ 
mulants. Similarly the cumulants determine the moments. 
In some cases theoretical treatments of problems in terms 
of cumulants are simpler than those of moments [Htta. 
The moment of a probability density function f{x) of 
a variable x is the mean value of x" and is mathematically 
defined by, 

+ CX3 

a„ = (x"-) = J x”/(x)(ix. (26) 

— OO 


r Cl = ai 


= a2 

- (ai)^ 




= ag 

— 3aia2 4- 

2(ai)^ 



= a4 

- 3 (a2)^ - 

- daiag -1- 

12(ai)"a2- 

■ 6(ai)'‘ 

= as 

— 5aia4 — 

10a2a3 -1- 

20a3 (ai)' 


+30 (3-2) 3-1 ■ 

-60 (ai)" 

32 + 24(31)^ 


= ag 

- Baias - 

153-23-4 + 

3O34 (31)^ — 

10 (ag)' 


-|-120aia2a3 - 120(ai)^a3 -|- 30(a2)^ 
-270(ai)2(a2)2 -h 360(ai)4a2 - 120(ai)6 
< C7 = a7 — 7aia6 — 21a2a5 -|- 42a3a5 — 35a3a4 

-|-210a4a2ai — 210a4a3 -|- 140aia§ -|- 210a3a2 

— 1260a3a2a3 -|- 840afa3 — 630a2ai -|- 2520afa2 
-2520afa2 -H 720a^ 

Cg = ag — 8aia7 — 28a2a6 -I- SGa^ag — SGagag 

-|-336a5a2ai — 336a5af — 35a4 -|- 560aia4a3 
-|-420a4a2 — 2520a4a2a^ — 1680afa4 -|- 560a|a2 

— 1680a^a3 — SOdOaga^ai -|- 13440a3a(^a2 
-6720afa3 - 630a4 -h 10080a2a3 - 25200afa2 
-H20160afa2 - 5040af 


(30) 
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' Mi = 0 

M 2 = a2 - (ai)^ = cr^ 

Ms = as - 3aia2 + 2 (ai)^ 

M 4 = a4 - daias + 6 (ai)^ a 2 - 3 (ai)"^ 

M 5 = as - 5a4ai - 10(ai)^a2 + lOaa (ai)^ + 4 (ai)® 
< Ale = ae - Gasai + 15(ai)2ji4 - 20a3 (ai)^ 

+15 (ai)^a2 - 5 (ai)® 

Mr = a? — Tagai + 21 ai^a 5 — 35a4af + 35a^a3 
— 21 afa 2 — 6a{ 

Ms = ag — 8a7ai + 28ai^a6 — SGagaf + 70a|a4 
—56afa3 + 28a3a2 — 7af 


(31) 

We can also write the cumulants in terms of central mo¬ 
ments: 

r Cl = ai 

C2=M2 

Cs = Ms 

C4 = M4-i{M2f 
C 5 = AI 5 — IOAI 2 AI 3 
^ Ce=Me- 15 M 4 M 2 - IQMl+iQMl 
Cr = Mr- 21 MsM 2 - ibM4Ms +210Ms Ml 
Cs=Ms- 28 MgM2 - 56 M 5 M 3 -35 Mj 
+ 42 OM 4 Ml + 5607^2 Ml - 6307W^ 


(32) 

which can be combined into a single recursive relationship 

n—1 

C„ = ^ C:i;-_\CmMn-m- (33) 

m—1 

General expressions for the connection between cumulants 
and moments may be found in m- A very convenient way 
to write the central moments and the cumulants in terms 
of determinants, 


C 


n 


and 


= (-1)”-' 


Ml 

1 

0 

0 

0 ... 

M 2 

Ml 

1 

0 

0 ... 

Ms 

M 2 

{T)Mi 

1 

0 ... 

M 4 

Ms 


q)Mi 

1 ... 

Ms 

M 4 

ibMs 


QMi ... 


(34) 


Cl 

-1 

0 

0 

0 ... 

C 2 

Cl 

-1 

0 

0 ... 

C 3 

(DC2 

()C3 

ri)C4 

Cl 

-1 

0 ... 

C 4 

C 5 

d)C3 

Cl -1... 
( 3)^2 Cl ... 


(35) 


where the determinants contain n rows and n columns 
and where (^) = standard binomial co¬ 

efficients. Then we can say that some important features 
of the system’s partition function can be deduced only by 
knowing all the moments. Each pf/i-order cumulant can 
be represented graphically as a connected cluster of p- 
points. If we write the moments in terms of cumulants by 


inverting the relationship (13011 or by expanding the deter¬ 
minant dsni), the pt/i-order moment is then obtained by 
summing all possible ways to distribute the p-points into 
small clusters(connected or disconnected). The contribu¬ 
tion of each way to the sum is given by the product of 
the connected cumulants that it represents. Due to the 
very important mathematical properties of the connected 
cumulants, it is often more convenient to work in terms 
of them. Henceforth and solely for simplicity, the word 
cumulant, implicitly means connected cumulant. 


3.2 Connected Cumulant Ratios Formalism 

In a symmetric distribution, every moment of odd order 
about the mean (if exists) is evidently equal to zero. Any 
similar moment which is not zero may thus be considered 
as a measure of the distribution’s asymmetry or skewness. 
The simplest of these measures is AI 3 , which is of the third 
dimension in units of the variable. In order to reduce this 
to zero dimension, and so construct an absolute measure, 
we divide by cr^. Reducing the fourth moment to zero di¬ 
mension in the same way as above we define the coefficient 
of excess (kurtosis) which is a measure of the fiattening 
degree of the distribution. In the literature, other expres¬ 
sions of skewness and kurtosis are used instead of what 
we have defined. Many other measures of skewness and 
kurtosis have been proposed (see for example Pearson in 

m)- 


3.2.1 General Definitions 


The Cumulants are considered as important quantities in 
physics but cumulant ratios are more important. Suggest¬ 
ing to review their definitions and deduce the most useful. 
Up to the present state of things, no formalism that would 
give the definitions of cumulant ratios in an unified way, 
exists. For this reason, it appeared to us instructive to try 
to standardize and unify the definition of the cumulant 
ratios in a clear and consistent way [15]. We start by the 
following definition : 




ncfn^r 


(36) 


which represents the generalized connected cumulant ratio 
between the cumulants {Cj} and the cumulants {C^} with 
positive exponents > 0 and (3j > 0} . From this 

definition we can distinguish four cases, namely : 


The Normalized Cumulant Ratios are obtained from (1551) 
when the following condition is fulfilled, 

^ a* X (i) = ^ Pj X (j). (37) 

i=i 4=1 
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The Unnormalized Cumulant Ratios are those ratios in 
which we have the contrary case, 

^ a* X (i) ^ ^ X (j). (38) 

i=l 3=1 

In this case we can distinguish two types of unnormal¬ 
ized cumulants: over-unnormalized cumulants in the case 
of ai X (i) < Pj X (j) and under-unnormalized cu- 

i=l 3 = 1 

mulants in the case of E ai X (i) > 

i=l 3=1 


The pt/i-Order Normalized Cumulant Ratios correspond 
to those in which only a pt/i-Order cumulant is suitably 
normalized : 

Vj G / Pj^p = 0 and Pp = I (39) 

thus 

(40) 

with 

n 

'^aiX{i)=p. (41) 

i=l 

The pt/i-Order Under-Normalized Cumulant Ratios which 
are the most useful ones. This time, we have a particular 
form of the latter case, in which the indices {i} are all less 
or equal to p : 



(42) 

i=l 


P 


Qfi X (z) = p, 

(43) 


i=l 


The numbers at are either integers or rational numbers. 
If we solve the last algebraic equation (H5|) , we obtain the 
values of {ai} for every definition. For example n = 4 : 

4 

ai X (i) = ai+ 2a2 + 3q; 3 -b Aa^ = 4. (44) 

i=l 


Table 1. Some p-Order Under-Normalized Cumulants 


pt/i-Order 


^<V 

2 

{ai = 2} 

^^1(01=2)1 ~2 

/C <2 — (y 

3 

{q 2 = 3/2} 

^^2 =3/2)1 ^ ^ 

4 

{02 = 2} 


5 

{q 2 = 5 / 2 } 

^^2=5/2)} ^ jj _ 


{02 = 1 , 03 = 1 } 

^{(.2 = i ).(.3 = i)> _ rj 
/V,<6 — 11 + 

6 

{02 = 3} 



{03 = 2} 

^^(.3=2)} ^ 


{02 = 1, 04 = 1} 


7 

{02 = 7/2} 

= Z71 


{02 = 1 , 05 = 1 } 

E(;-^'(-)>^.2 


{02 = 2, 03 = 1} 

4(;-)’(--)>^.3 


{03 = 1 , O4 = 1 } 

;.{(.3 = l)-(<»;=l)l 
/C<7 * = T ] 4 , 

8 

{02 = 4} 

= u \ 


{02 = 1, 06 = 1} 

;.{(.2 = i)’(<»6 = i)> 

/C<g ^ ^ = UJ 2 


{02 = 2,04 = 1} 

/V<'g — CLI3 


{03 = 1,05 = 1} 

/C<g =0)4 


{02 = 1, 03 = 2} 

;.{(.2 = i)’(. 3 = 2 )} 

/C<g ^ =0^5 


{04 = 2} 

;.f(.;=2)} _ 

/v-^^g — L(J6 


with 


When solving this equation in the set of natural numbers 

N we hnd only five possibilities: E ^ “ E ^ ~ P' 

{(ai, a2, as, a4) = (0, 2, 0,0); (1,0,1,0); (0,0, 0,1); 

(2,1,0,0);(4,0,0,0)}.(45) 


From the relations (pTH]) and (HU) we derive the relation¬ 
ship which combines two different definitions of pt/i-order 

under-normalized cumulant and 

, which is given by. 




(46) 


From the relation we see that the number of possible 

definitions of increases with the order p. How¬ 

ever, we shall not consider all definitions, but we focus 
only on those mostly used. Generically, the structures of 
all cumulants are related to each other and the behavior 
including the magnitudes can be deduced from the pre¬ 
ceding. 























M. Ladrem et al.: Finite Volume Cumulant Expansion in QCD-Colorless Plasma 


3.2.2 The First Order Under-Normalized Cumulant Ratio: 
Normalized Mean Value 

Because the first cumulant is the mean value of a: :Ci = 
ai = (x) then the first order under-normalized cumulant 

ratio is = 1 


3.2.3 The Second Order Under-Normalized Cumulant Ratio 
: Normalized Variance 


The second order under-normalized cumulant ratio may 
be referred to the normalized variance and defined as, 


K. 




<2 


{xY 


9i 

c? 


M2 

Cl 


(xy 


- 1. 


(48) 


3.2.4 The Third Order Under-Normalized Cumulant Ratio: 
Skewness 


The third order under-normalized cumulant ratio is a mea¬ 
sure of symmetry, or more precisely, the lack of symme¬ 
try. A distribution, or data set, is symmetric if it looks 
the same to the left and the right of the center point. The 
3’^'^ cumulant for a normal distribution is zero, and any 
symmetric distribution will have a third central moment, 
if defined, near zero. Then the 3'''^ under-normalized cu¬ 
mulant ratio is called the Skewness S and is defined as. 


K 


{U4/2)} 


<3 


= A = 


(C2) 


3/2 


Ms 

mT' 


(49) 


A distribution is skewed to the left (the tail of the distri¬ 
bution is heavier on the left) will have a negative skewness. 
A distribution that is skewed to the right (the tail of the 
distribution is heavier on the right) will have a positive 
skewness. 


3.2.5 The Fourth Order Under-Normalized Cumulant Ratio: 
Kurtosis 


The fourth order under-normalized cumulant ratio is a 
measure of whether the distribution is peaked or flat rel¬ 
atively to a normal distribution. Since it is the expec¬ 
tation value to the fourth power, the fourth central mo¬ 
ment, where defined, is always positive. Because the fourth 
cumulant of a normal distribution is 3 ct^, then the most 
commonly definition of the fourth order under-normalized 
cumulant ratio called Kurtosis k is 


K. 




<4 


= K = 


C4 

(C 2 )" 


Mi 

Ml 


3, 


(50) 


so that the standard normal distribution has a kurtosis of 
zero. Positive kurtosis indicates a ’’peaked” distribution 
and negative kurtosis indicates a ’’flat” distribution. Fol¬ 
lowing the classical interpretation, kurtosis measures both 
the ’’peakedness” of the distribution and the heaviness of 


its tail m- In addition to this. Binder was the first who 
has proposed and studied the fourth cumulant as it was 
defined in |21j using the moments of the energy probability 
distribution: 


^4 = 1- 


sg 

3 al' 


(51) 


This was introduced as a quantity whose behavior could 
determine the order of the phase transition. If we replace 
the moments by the central moments, we get another com¬ 
pletely different physical quantity, which is related the 
kurtosis as. 


Mi 
3Ml ~ 3 


(52) 


and can be easily derived from our general definition of 
connected cumulant ratios dsni)- This new cumulant, as we 
have mentioned before, is called connected Binder cumu¬ 
lant or conventional Binder cumulant. However, to avoid 
confusion in the appellations we simply keep the name of 
Binder cumulant for the first quantity. Historically, this 
new cumulant was first introduced and studied by Binder 
in 1984 [5D]. Seven years later, this new cumulant was 
reconsidered in an independant and important paper by 
Lee and Kosterlitz in the context of a different model 
|22) . The difference between the two Binder cumulants 
attracted little attention in its early years. But, in 1993, 
Janke has illuminated the most important difference in 
a comparative and fruitful study between the two cumu¬ 
lants [23]. The great significance of the connected Binder 
cumulant relative to the Binder cumulant is summed up 
in the following points: (1) the thermal behaviors of two 
Binder cumulants are very different, particularly in the 
transition region, (2) the connected Binder cumulant has 
a richer structure than the Binder cumulant, (3) the con¬ 
nected Binder cumulant is more efficient in locating the 
true finite volume transition point than the Binder cumu¬ 
lant. This cumulant is a finite size scaling function |2ni 
and is widely used to indicate the order of 
the transition in a finite volume. In ordered systems, a 
good parameter to locate phase transitions is exactly this 
connected Binder cumulant which is the kurtosis of the 
order-parameter probability distribution. The uniqueness 
of the ground state in that case is enough to guarantee 
that the Binder cumulant takes the universal value at zero 
temperature for any finite volume. 


3.2.6 The Fifth Order Under-Normalized Cumulant Ratios: 
Pentosis 

The fifth order under-normalized cumulant ratio, which is 
called Pentosis, can be defined in two ways. The first one 
is given by. 


K 




= n+= 


Mr 


C2C3 M2M3 

and the second definition is given by, 

C5 Ms 


- 10 (53) 


^^(.2=5/2)} 


(C2) 


5/2 


M 


5/2 


- lOA. 


(54) 
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The two forms of pentosis are of course related by, 


3.2.10 Three Unnormalized Cumulant Ratios 


77_ = Sn+. (55) 

a relation that we can deduce from the general relationship 

(HU). 


We are also interested in studying different unnormalized 
combinations of the cumulants. Their importance was re¬ 
vealed and emphasized in several recent works Esiisniisii 
132] . The first combination contains the variance cr^, kur- 
tosis K and skewness S and is defined as, 


3.2.7 The Sixth Order Under-Normalized Cumulant Ratios: 
Hexosis 

The sixth order under-normalized cumulant ratio is, anal¬ 
ogously to Pentosis and Kurtosis, coined Hexosis. It can 
be dehned in one of the following ways mm- 


^ ^ _ C^/^C4 _ _ -Ml {M4 - 3M1) 

S Cs {( 3 )} M 3 

(62) 

The second one contains only the variance tr^ and skew¬ 
ness S and is given by. 


K 


{(..U)} 


<6 


= 77i = 


Mq - 157W4M2-107U1§ 


(Cz)^ 


Ml 


-h30, 

(56) 


U = 


1 


1 I /2 


r2S 


= /c 


{(1/=)} _ Ml 

fG)} 


M 3 


(63) 


K 


{{<, 3 = 2 )} _ T, _ Ce Me — X^M^Mi 'iQMl and kurtosis k and is given by. 

Ml 

(57) 


However,the third combination contains the variance 

Ma - 3Ml 


M = a^K = ^ = /C 


K 


{(..GG.UtOi 


<6 


= % = 


Me - IQMl-lbMl 


C4C2 


MA-3{M2y 


-15 


C 2 


fG)} 

fG)} 


M 3 


(64) 


M 3 


(58) 

It is easy to show that the three definitions of hexosis are 
related each other by the relations. 


3.3 Finite size cumulant expansion of the hadronic 
probability density function p(h) as function of 

L^n{q,T, V) 


He = (59) 

which can be deduced from the general relationship (H51) . 


3.2.8 The Seventh Order Under-Normalized Cumulant 
Ratios: Heptosis 

With the same spirit and by analogy to pentosis, kur¬ 
tosis and hexosis, we can term the seventh order under¬ 
normalized cumulant ratio as Heptosis m rj. One of the 
possible definition of heptosis is given by. 


K 


{(./=2).(.3Gl)} 


<7 


= % = 


M, 


- 21 - 


C^Ca 

Me 


MiMa M 2 M 3 


- 354^1 + 210, (60) 

Ml 


3.2.9 The Eighth Order Under-Normalized Cumulant 
Ratios: Octosis 


Concerning the eighth order under-normalized cumulant 
ratio, which can be termed Octosis m and take one of 
the eight definitions from table o. 


K 




<8 


-56 


Me 


— ue — 


Ml 


C 2 CI 


•^8 _ 00 Ale 

MIM 2 Ml 

.Ml 


M 3 M 2 


-35-^^ +420-630:^^+560. 

MIM2 Ml Ml 


Using our hadronic probability density function p(h), we 
derive the general expression of the mean value (h") as 
a function of Lmn{q,T,V) [33]- Afterwards, one can ex¬ 
press the different cumulants C„(T, U) in terms of these 
Lmn{q,T,V) using (l2^ and (l35ll . Keeping in mind that 
these double integrals Lmn{q^ T, V) are state functions de¬ 
pending on the temperature T, on the volume V and on 
the state variable q. One can hide their dependence on 
(T, V) just to avoid overloading relationships. After some 
algebra, we get the result. 


(h”) (T V) = (fc)fc!Ao^fc-|_i (0) 

Ao.i (1) — To.i (0) 

(65) 

Using this general expression of the mean value and from 
dSOl) we derive the six first cumulants (see appendix), 

rCi(r, u) = (h) 

C2(T, U) = (h2)-(h)' 

C3(T, U) = (h3)-3(h)(h2) + 2(h)3 
C 4 (r, V) = (h^) - 3 (h^)" - 4 (h) (h3) +12 (h)' (h2) 
-6(h)" 

Ce{T, V) = (h5) - 5 (h) (h4) - 10 (h2) (h3) 

< +20(h3)(h)V30(h2)^(h) 

-60 (h)^(h2) +24(h)5 
C6(T, U) = (h6) - 6 (h) (h5) - 15 (h2) (h4) 

+30 (h^) (h)^ - 10 + 120 (h) (h^) (h^) 

-120 (h)^ (h3) + 30 (h2)^ - 270(h)2(h2)2 
+360(h)4(h2) - 120(h)6 


(61) 


(66) 
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Afterwards we derive the final expression of both p-order 
under-normalized and unnormalized cumulants under con¬ 
sideration. The first cumulant is none other than the order 
parameter (h) {T,V) and is given by (1^^ . The variance 
cr^ (T, V) is given by, 


<7^ (T,V)= [(h2)-(h)2] . 

The skewness S {T, V), is given by, 

S (T V) = = [(h3)-3(h)(h^) + 2(h) 

’ [(h2) - (h)2]3/2 

and the kurtosis n (T, V)is given by. 


(67) 


( 68 ) 


.(T,V)=(("-lh))^)-3 = 


[(h4) - 4(h)(h3) - 6(h)4 + 12(h)2(h2) - 3(h2)2 

[(h^) - (h)^]" 


(69) 



Fig. 1. Behavior of Different Cumulants C„=i, 2 , 3 , 4 , 5 , 6 (T, V) vs 
Temperature for Volume = 1000/m®. 


And finally the pentosis 11+ (T, V), which is given by, 


n+ (T, V) = Nl/Dl 

N1 = [(h5) - 5(h4)(h) + 20(h3)(h)2 - 60(h2)(h)3 
< -10(h2)(h3)-h30(h2)2(h)-t24(h)^] . 

D1 = [(h2)(h3) - 3(h)(h2)2 + 5(h2)(h)3- 

(h)2(h3)-2(h)5] 

(70) 

And finally the hexosis V), which is given by, 

-Hi (T, V) = N2/D2 

N2 = [(h6) - 6(h5)(h) - 15(h2)(h4) -h 30(h)2(h4) 

, -10(h3)2 -h 120 (h) (h2) (h®) - 120 (h)^ (h®) 

-f30 (h^)^ - 270(h)2(h2)2 -p 360(h)4 (h^) - 120(h)6 
^ D2 = [(h2)3 + 3(h)4(h2) - 3(h2)2(h)2 - {hf] 



T[MeV] 


(71) Fig. 2. Hadronic Probability 

The expressions of 77_ {T,V) and T-L 2 ,d,{T,V) can be de- p(h,T, V = 1000/m® 
rived easily from those of (1551) and (1551) using (1551) . Let us ues of h = 0.1,0.5,0.9 
now go to the unnormalized cumulants as defined in ((621631641) . 

Final expressions of 0,U,M are : 


Distribution Function 
vs Temperature for different val- 


O (T, V) = 


(T, V) AC (T, V) 


s(r,v) 


[(h4)-4(h)(h3)-6(h)^ 


-bl2(h)2(h2)-3(h2)2] 

[(h2)-(h)2]^/^ 


U{T,V) = 


and 


[(h3)-3(h)(h2)+2(h)3] 
[(h^) - (h)2] 


(72) 


1/2 


M{T,V) = a^ (T, V)k{T, V) = 

[(h4) - 4(h)(h3) - 6(h)4 + 12(h)2(h2) - 3(h2)2] 

[(h^) - m 


(74) 


We will see after studying these new thermodynamic func¬ 
tions, that their FSS analysis will allow to identify the 
transition region, to define judiciously the finite volume 
transition point and analyse its behavior when approach¬ 
ing the thermodynamic limit. 


4 Finite Size Cumulant Expansion 
and Discussion 


Results 


a2 (T, V) S (T, V) [(h3) - 3(h)(h2) + 2(h)3] 

(73) 


Firstly, one may notice a clear sensitivity, of the whole 
quantities studied in this work, to finite volume of the 
system. Exactly as in the case of the results obtained in 
our previous work mu. the 4FSE cited above are ob¬ 
served. The variation of the different cumulants and cu¬ 
mulant ratios versus temperature are illustrated in Figs. 
(1-13) respectively, for various finite sizes. They show in¬ 
teresting features. It can be clearly seen that the differ- 
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Fig. 3. 3-Dim plot of the Order Parameter < h(r, V) > vs 
Temperature and Volume. 

Fig. 4. 3-Dim plot of Variance (T, V) vs Temperature and 
Volume. 

Fig. 5. 3-Dim plot of Skewness S (T, V) vs Temperature and 
Volume. 

Fig. 6. 3-Dim plot of Kurtosis k {T, V) vs Temperature and 
Volume 

ent finite peaks appearing in the different quantities have 
width 5T [V) becoming small when approaching the ther¬ 
modynamic limit. This result is expected, since the order 
parameter looks like a step function when the volume V 
goes to infinity, as it is well known. The rounding of the 
cumulants behavior is a consequence of the finite size ef¬ 
fects of the bulk singularity. We notice in all curves, the 
emergence of a transition region, roughly bounded by two 
particular points, which narrows as the volume increases. 
In this region, all thermodynamical quantities present an 
oscillatory behavior which becomes faster when approach¬ 
ing the thermodynamic limit. Our previous works 
have shown that both (h) and exhibit a finite sharp 
discontinuity, which is related to the latent heat of the 
deconfinement phase transition, at bulk transition tem¬ 
perature To (oo) = = 104.34796MeV, reflect¬ 

ing the first order character of the phase transition. It 
is well known that the latent heat is the amount of en¬ 
ergy density necessary to convert one phase into the other 
at the transition point. In our case, the latent heat can 
be calculated:T// (oo) = 45. This finite discontinuity can 
be mathematically described by a step function, which 
transforms to a d-function in xt and ct- When the vol¬ 
ume decreases, all quantities vary continuously such that 
the finite sharp jump is rounded off and the d-peaks are 
smeared out into finite peaks over a range of temperature 
6 T {V). Physically, we can interpret these 4FSE as due to 
the finite probability of presence of the CQGP phase below 
the transition point and of the hadron phase above it, in¬ 
duced by the considerable thermodynamical fluctuations. 
In Fig. (I), we show the plot of the first six cumulants 


as functions of temperature at fixed volume :1000/m^. A 
multiple peaks structure can be observed on these curves, 
except in the case of the first cumulant C'i(T). For each 
additional order, a new hump (peak)is introduced. These 
peaks are broadened, smaller is the volume. Also, we no¬ 
tice that the inflection point in the first cumulant C'i(T) 
becomes a maximum point for the second order cumulant 
(72(T), a zero point in the third cumulant C^lT) and so 
on. The number of times that a given cumulant changes 
its sign, is directly related to the order of the cumulant. 
The sign change for the cumulants starts at the third one. 
It happens twice in the fourth, thrice in the fifth and four 
times in the sixth order cumulants. The common feature 
is that the higher the order of the cumulant, the higher 
frequency of the fluctuation pattern is. Also, we notice 
that all cumulants have the same vanishing value at low 
or high temperatures. In the middle region, which in prin¬ 
ciple is considered as the transition region, the value of 
the cumulants presents an oscillatory behavior due to the 
thermodynamical fluctuations during the phase transition. 
When we carefully analyse the behavior of the hpdf for 
different values of h = 0.1,0.5 and 0.9 on Fig.2, we note 
that in the case of h = 0.5 the hpdf looks like very sym¬ 
metric and for these reasons we expect the skewness to 
be zero. The hpdf distribution is skewed right before the 
transition h = 0.1 and becomes skewed left after the oc¬ 
currence of the phase transition h = 0.9. The peaks of 
the hpdf are more pronounced when we go from a pure 
CQGP phase to a pure hadronic phase passing throught 
the mixed phase.This feature is simply due to the fact that 
our hpdf is directly connected to the density of states in 
each phase. 

Let us now see what the plots of the normalized cumu¬ 
lants in Fig. (3-13) express? The general behavior and the 
structure of the peaks are much different. However the 
broadening effect of the transition region with decreasing 
volume is also observed. The plots of skewness, kurtosis 
and pentosis, show a double peaks structure, a big peak 
and a little one. These two peaks correspond to the two 
states before and after the phase transition. When the two 
peaks have the same sign, there are two vanishing points 
limiting the transition region and containing a small ex¬ 
tremum which is nothing other than the transition point. 
This behavior is due to the fact that kurtosis is closely 
connected the second derivative of the thermal suscepti¬ 
bility. Otherwise there is only one vanishing point which 
is the transition point. The only difference between the 
three curves lies on the fact that the small peak becomes 
less pronounced with increasing order of the cumulant. 
For this reason, the latter does not appear practically on 
the curves. In the transition region the symmetric peak of 
p{h = 0.5, T) becomes very small by making the kurtosis 
negative and small. The kurtosis manifests a very differ¬ 
ent behavior in both sides of the transition region when 
approaching the thermodynamic limit which is due to the 
high asymmetry of the variance, as displayed clearly on 
the 3-Dim plot in Fig. (4). The variance decreases more 
sharply in the hadronic phase than in the GQGP phase. 
When looking more closely at all the 3-dimensional plots, 
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we can clearly see that some particular points exhibit a 
typical behavior that can be described by the finite size 
scaling law, which is consistent with what has been ob¬ 
tained previously [5]. For example, the maximum of the 
variance, sketches the finite size scaling behavior described 
by : r((Tmax) — 70 ( 00 ) oc V~^. Concerning the plots of the 
three hexosis, namely 'Hi, 2 , 3 , a same global behavior out 
of the transition region and a different oscillatory behavior 
in it. The local maximum point in "Hi becomes a singular¬ 
ity point in "^2 and a local minimum in 'Ha. Moreover, the 
obvious change in the sign, observed in our results, is in 
agreement with the results obtained by other models m- 
Finally the plots in Fig.11,12,13 represent the variations of 
the three unnormalized cumulant ratios O {T, V ), U (T, V) 
and Af (T, V) as a function of temperature and volume. 
Their behaviors are very different compared to the plots of 
the normalized ratios. The plots of J\f (T, V) show a clear 
and rapid oscillatory behavior with two maxima and one 
minimum in the transition region which gradually narrows 
as the volume increases. On the other side we can clearly 
see the emergence of particular singular behavior on the 
plots of O (T, V) and U (T, V) at certain values of tem¬ 
perature. The same divergence is observed on the plot of 
the pentosis n+(T, V), exactly in the valley region be¬ 
tween the two maximums (Fig.9). It is interesting to note 
the behaviour of O (T, V) which is practically zero in the 
two phases and is singular at the finite volume transition 
point, with a small local minimum before the transition 
and small local maximum after the transition.The location 
of the finite volume transition point is clear and simple, 
its shifting is obvious. The same observations are valid 
for U {T,V). Using FSS analysis, we will see below that 
these points will be identified as the finite volume tran¬ 
sition points. We summarize by saying that O (T, V) and 
U (T, V) tend to zero rapidly everywhere, except in the 
transition region and at the finite volume transition point 
where they diverge. This, is due to the zero of skewness in 
the transition point. These two cumulant ratios can there¬ 
fore serve as two good indicators of the location of the 
finite volume transition point. They will be of great use 
in the analysis of experimental data of URHIC where the 
context of initial conditions just before the phase transi¬ 
tion are unknown. We can see again from the figures that 
change their values sharply from negatives to positives 
and oscillate greatly with temperature near the transition 
point. These qualitative features; ie, sign change and oscil¬ 
lating structure, are consistent with effective models [5^ . 


5 New Method of Localization of the Finite 
Volume Transition Point 

5.1 Natural Method 

It is important to have a precise knowledge of the region 
around the transition point since many quantities of phys¬ 
ical interest are just defined in its vicinity. It therefore 
seems very important to find the definition of a finite- 
volume transition point which involves less corrections. 



Fig. 7. Pentosis II_ (T, V) vs Temperature for different Vol- 
umes(-|- zoom of the region close to zero). 



Fig. 8. Pentosis 11+ (T, V) vs Temperature for different Vol¬ 
umes (-1- zoom of the region close to zero). 



Fig. 9. Pentosis 11+ (T, V) vs Temperature for Volume = 
100/mT 
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Fig. 10. Hexosis 'Hi, 2,3 {T,V) vs Temperature for Volume — Fig. 11. 0{T,V) vs Temperature for different Volumes (+ 
100/m® (+ zoom of the region close to zero). zoom of the region close to zero). 


Let us first remind the logical and natural way to define 
the finite volume transition point by saying that is the 
point where we have equal probabilities between hadronic 
phase and CQGP phase: (h (To"(P))) = l-(h (To^(V))). 
This means that the value of the order parameter is given 
by (h {T^(V))) = 1/2. We know that in thermodynamic 
limit the order parameter manifests a finite discontinu¬ 
ity which can be easily described by a step function (El- 
Therefore, the specific heat ct {T, V) and the thermal sus¬ 
ceptibility XT {T, V) show i5-function singularities at the 
transition point : 


lim I I oc 5{T — To(oo)) (75) 

(y)^oo \ XT (P, P) j 

In hnite volume, these ^-singularities become rounded peaks. 
Therefore xt (P, V) and ct (T, V) reach a local extremum 
value at certain temperature (V) which is defined as 
the temperature of the hnite volume transition point : 

{ xt(t'v) =»“:}^(™) 

Finally, we can assert without any problem that the hnite 
volume transition point (and its temperature T^{V)) is 
logically the point where the following equations are sat- 
ished : 


dxT(T,V) 

dT 


(h(To^(I/))) = l/2 
= 0 and 

TN (y) CO 


To^iV) 


= 0 


(77) 



Fig. 12. U (T, V) vs Temperature for different Volumes. 



Fig. 13. M {T, V) vs Temperature for different Volumes. 


From this we see that the hnite volume transition point 
is associated to the appearance of an inhection point in 
(h(T, V)): (h {T^{V))) becoming a local extremum point 
in both XT {T, V) and ct (T, V). According to this method, 
we extract the different temperatures (V) of the tran¬ 
sition points and collect them in table ©■ 
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Table 2. Natural Transition Points Temperatures 


Volume V fm^ 

(V) [MeV] 

100 

110.68007 ±0.00001 

200 

108.02068 ± 0.00001 

300 

107.00271 ± 0.00001 

400 

106.44758 ± 0.00001 

500 

106.09471 ± 0.00001 

700 

105.66892 ±0.00001 

900 

105.41823 ±0.00001 

1000 

105.32709 ± 0.00001 

2000 

104.88739 ± 0.00001 

5000 

104.59347 ± 0.00001 

10000 

104.48254 ± 0.00001 

oo 

104.34796 



Fig. 14. Correlation scatter plot between ro(F){Qn,} + 
n{2MeV) and Tq {V) for different volumes {Qo = o-'^ax^Qi = 
— lamin'; Q3 = (17+)oo,g4 = (1 T_)o,Q 5 = OocQe = 
U oo ^ Qy — 

5.2 Cumulant Method : Particular Points and 
Correlations 

In this section, we will try to propose a new method for 
locating the finite volume transition point using the whole 
cumulants studied in this work. We shall show how this 
finite volume transition, clearly manifests itself as a partic¬ 
ular point in each cumulant. Our strategy, we use, consists 
of finding a judicious point where the temperature T^iy)^ 
seemingly tends to the bulk To(oo) with increasing volume 
and must be highly correlated with {V) : 

lim ro(V) = ro(oo). (78) 

(V)-)-oo 

The definition of ro(fo) is not arbitrary but very difficult 
analytically and differs according to the quantity being 
considered. After a careful analysis of the normalized cu¬ 
mulants plots cr2(r, V), E{T, V), k{T, V), n±{T, V), 
ni,2,3{T,V),0{T,V),U{T,V) and Af{T,V), we find that 
the only points which can be considered in one way or 
another as very particular are : the local extrema points 
(local maximum and local minimum), the vanishing points 


(zeros), the inflection points and the singular points. These 
points are called the Particular Points. Indeed, we have in¬ 
vestigated the behaviour of these particular points. Firstly, 
for each quantity and for each particular point, we extract 
the temperature values {To(fo)} at different volumes and 
put them in the first set. Secondly we put the tempera¬ 
ture values {To^(fo)} given in table in the second set. 
To probe more precisely the location of the finite volume 
transition point, a useful tool is the scatter plot, in which 
the temperatures of the first set are plotted against the 
temperatures of the second set. What we are asking here 
is whether or not the variations in the first set of Tq{V) 
are correlated or not with the variations in the second set 
of {V). We have analyzed several particular points and 
only good candidates are considered in this work with de¬ 
tails. If a particular point is considered as a good finite 
volume transition point, one would expect that its scatter 
plot satisfies the following three criteria: 

(1) The fit should be linear. 

(2) The slope of the fit should equal unity and its vertical 
intercept should equal zero. 

(3) The fit should have high linear correlation with a very 
good correlation factor and a very good probability test. 
If we consider the temperature {To{V)} to be dependent 
variable, then we want to know if the scatter plot can be 
described by a linear function of the form, 

ToiV) = Xy^{V) + ,y. (79) 

Because we are discussing the relationship between the 
variables {To(I7)} and {Tj^(fo)}, we can also consider 
{Tj^(t7)} as a function of {To(fo)} and ask if the data 
follow the same linear behavior, 

T^{V) = X'ToiV) + i2'. (80) 

The values of the coefficients X' and z/' in dHOD will be 
different from the values of the coefficients A and ly in 
equation GSD, but they are related if the two temperatures 
{To{V)} and {TJ^ (1^)} are correlated. If we consider solely 
the value of A (or A'), it doe not provide us a good measure 
of the degree of the correlation. From (1751) and (|5(I)) . and 
in the case of a total correlation, we can show that 

If there is no correlation, the two parameters A and A^ 
are lower than unity, even approaching zero value. We 
therefore can use the product AA' as a measure of the 
correlation between the two sets of temperatures {To{V)} 
and (1^)}- By definition the correlation factor is given 
by p = VXX'. The value of g ranges from 0 , when the 
data are totally uncorrelated, to 1, when there is total 
correlation. The correlation factor, alone, is not sufficient 
to indicate the quality or the goodness of the linear fit. 
An additional calculation of probability is necessary for 
more precision. This probability distribution enables us 
to go beyond the simple fit, and to compute a probability 
associated with it. In the case of our situation, a commonly 
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Table 3. Correlation factor values obtained from linear fitting 


N.Cumulant 

Transition Point 

A 

A^ 


<.AToiV)) 

0.98812 

1.01202 

r(r, F) 

So{To{V)) 

0.98798 

1.01216 

k{T,V) 

«:min(To(F)) 

0.98700 

1.01317 

n+{T,v) 

in+)^ (Toiv)) 

0.98788 

1.01226 

n.{T,v) 

(77), (To(F)) 

0.98753 

1.01262 

0{T, V) 

Goc(ro(F)) 

0.98787 

1.01227 

U{T, V) 

Uao{To{V)) 

0.98787 

1.01227 

J\f(T, V) 

MminiMV)) 

0.98753 

1.01262 


used probability distribution for g is given by [36l|37], 

(o C) — —- p2')(C-2)/2 /•g2') 

I ^ r[{C)/ 2 ] ^ 

where ( = N — 2 is the number of degrees of freedom 
for a sample of N data points, and r{x) is the standard 
Gamma function.lt gives the probability that any sam¬ 
ple of uncorrelated data would yield to a linear behavior 
described by a correlation factor equal to g. If this prob¬ 
ability is small, then the sample of data points can be 
considered as highly correlated variables. More generally, 
this type of calculation is often referred to as goodness of 
ht test |38] . Another significant and useful quantity which 
can be calculated from the distribution (15^ is given by, 

1 

Vc{g,N) = 2jv,{xX)dx. (83) 

hi 

This Pc{g,N) represents the integral probability that a 
sample of N uncorrelated data points would yield a linear 
correlation factor larger or equal than the calculated value 
of |p|. This would mean that a small value of PciQ, N) is 
equivalent to a high probability that the two sets of vari¬ 
ables are linearly correlated. The fitting results obtained 
from the correlations study showed on Fig.(14), are sum¬ 
marized in table Q . In order to avoid overlapping between 
fitting curves and to allow a clear representation on the 
same graph, we have added a shift of 2 MeV between each 
two consecutive curves. It can be perceived from the scat¬ 
ter plots Fig. 14 that the points are closely scattered about 
an underlying straight line, refelecting a strong linear re¬ 
lationship between the two sets of data and the numerical 
values of the slopes are close to unity as expected. Also, 
we tried the fitting procedure with a fixed intercept v = 0 
and we got better results, the value of the slope better 
than 0.999. From the values of both A and A' in the table 
3, pratically a same value of the correlation factor p, which 
equal to 0.99999, is obtained. Therefore the evaluation of 
the two probabilities gives the following results: 

P^ig = 0.99999, ( = 7) = 1.82209 x IQ-^^ 

Pc(g = 0.99999, N=9) = 1.04119 x IQ-^^ ’ 

The extreme smallness of Pc{g,N) < 1.178 x 10“^® in¬ 
dicates that it is extremely improbable that the variables 


under consideration are linearly uncorrelated. Thus the 
probability is very high that the variables are correlated 
and the linear fit is justified. The fact that such fittings 
yield results that are consistent with each other is an im¬ 
portant consistency check on the accuracy of the calcula¬ 
tions and gives an idea of the FSE for the values of the 
temperature of finite volume transition point . We would 
like to note that the numerical values of temperature ob¬ 
tained by the cumulant method {To(k^)} of the various 
transition points, are comparable with an accuracy less 
than 2%, with the temperatures (1^)} extracted using 
conventional procedures. Therefore the selected points are 
indeed the true finite volume transition points, namely: 

(1) the local maximum point in the variance cr^(T', V) and 
in the first hexosis 'Hi{T,V): 

(2) the zero point in the skewness S{T, V) and in the pen- 
tosis n-(T,V): A7o,i7_^o, 

(3) the local minimum point in the kurtosis k{T,V), in 
Af{T, V) and in the third hexosis TLsiT, V): /Cmin, ^ 3 .mi™, 
Amin 

(4) and the singularity point in the pentosis n^{T, V), in 

i({T, V), in 0(T, V) and in the second hexosis V): 

n+,coi Uoo. O OO j ’^2,00- 

The temperature at which skewness vanishes is expected 
to represent the transition temperature, and tends appar¬ 
ently to Tq{oo) with increasing volume, while the temper¬ 
ature gap between the two extrema is expected to give the 
width of the transition region. 

We got an unexpected and important result. It concerns 
the behavior of the connected Binder cumulant. Indeed 
from the relation(34) the whole discussion about the kur¬ 
tosis can be translated to the connected Binder cumulant. 
Therefore, the connected Binder cumulant B%{T,V) has 
two minima and a little maximum between them as ex¬ 
pected from the behavior of the kurtosis k{T, V). The po¬ 
sition of two minima should not have a good correlation 
factor, however, the little maximum will be the good fi¬ 
nite volume transition point. This would be in striking 
contrast to conventional result obtained by Binder m- 
The apparent discrepancy is completely due to the differ¬ 
ence in the defintion of the Binder cumulant B^iT^ V) and 
the connected Binder cumulant B%{T^V). The local min¬ 
imum point in the Binder cumulant is not the true finite 
volume transition point because it has not the good corre¬ 
lation factor (A = 1.39). But it should approach the bulk 
transition temperature as V becomes large, which means 
that it is just a particular point. We have therefore shown 
that the cumulants are more interesting than the moments 
and the connected Binder cumulant is more efficient in 
locating the true hnite volume transition point than the 
Binder cumulant. The same results are obtained in many 
papers and the obtained thermal behaviors are 

in complete agreement with ours. We know that all the 
particular points as they have been defined in our paper 
converge towards the unique singularity in the thermody¬ 
namic limit. Once the true finite volume transition point 
has been identified from the particular points, its signal 
is not necessarily the highest, and even, may be in some 
cases, is hard to detect. The main property of the particu- 
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lar points in finite volume is that they are correlated with 
the true finite volume transition point. Another important 
property relates to the possibility of using them to define a 
transition region. It has been claimed in that the shift be¬ 
tween the minimum of the Binder cumulant and the max¬ 
imum in its susceptibility in the case of a first order phase 
transition, is due to the absence of the phase coexistence 
phenomena in the double Gaussian model and of the sur¬ 
face corrections [^[^ . In our case, despite the intake into 
account of the phase coexistence within the Colorless-MIT 
bag model, the shift between the minimum of the Binder 
cumulant and the true finite volume transition point still 
exists but its magnitude is different. The magnitude of 
this shift is reflected in the numerical values of the corre¬ 
lation parameters (A, v) which differ from the ideal values 
(A = l,z^ = 0) in the case of a total correlation. Indeed, 
when we try to extract roughly the numerical values of 
A parameter from the results obtained in 
we find different values [A = 1.55,1.47,1.57,1.89] respec¬ 
tively, which are not close to unity. This is certainly due to 
the fact that our Colorless-MIT bag model is very different 
from the double gaussian model used by Binder to study 
the finite size effects in the first order phase transition EH. 
Presumably the shift of the minimum of B 4 {T, V) from the 
true finite volume transition point (V) depends on the 
detailed form of the partition function of the system un¬ 
der consideration as quoted in |22j . ie, it is somewhere 
model-dependent. 

6 Conclusion 

In order to identify and locate the finite volume transi¬ 
tion point more accurately, we have studied in details the 
finite volume cumulant expansion of the order parameter 
and have shown how greatly this can be used to provide 
a clear definition of the finite volume transition point in 
the context of the thermal deconfinement phase transi¬ 
tion to a CQGP. Starting from the hadronic probability 
density function and using the Lmn-method, a finite size 
cumulant expansion of the order parameter is carried out. 
The first six cumulants, their under-normalized ratios and 
also some combinations of them, are then calculated and 
analyzed as a function of temperature at different vol¬ 
umes. To be more consistent and coherent in our defi¬ 
nitions of cumulant ratios, a new reformulation of these 
cumulant ratios is proposed. It has been put into evidence 
that all cumulants and their ratios showed deviations from 
their asymptotic values (low and high temperature values), 
which increase with the cumulant order. This behavior is 
essential to discriminate the phase transition by measur¬ 
ing the fluctuations. We have noticed that both cumulants 
of higher order and their ratios, associated to the thermo¬ 
dynamical fluctuations of the order parameter, in QCD 
behave in a particular enough way revealing pronounced 
oscillations in the transition region. The sign structure 
and the oscillatory behavior of these in the vicinity of the 
deconfinement phase transition point might be a sensitive 
probe and may allow to elucidate their relation to the 
QCD phase transition point. In the context of our model, 


we have shown that the finite volume transition point is 
always associated to the appearance of a particular point 
in whole cumulants under consideration. A detailed FSS 
analysis of the results has allowed us to locate the finite 
volume transition points and extract accurate values of 
their temperatures To{V). We have tested the validity of 
our results by performing linear correlations between the 
set oiTo{V) and the known results obtained with the nat¬ 
ural definition {V) providing very good correlation fac¬ 
tors. In addition to natural definition of the finite volume 
transition point as the extrema of thermal susceptibility, 
XT and specific heat ct, we have shown that the true finite 
volume transition point manifests itself as a different par¬ 
ticular point according to the quantity considered, namely 
as, 

(1) a local maximum point in the variance cr^(T, V) and 
in the first hexosis 'Hi{T,V): Cmax! 

(2) a zero point in the skewness S{T, V) and in the pen- 
tosis n-{T,V): A'o,i7__o, 

(3) a local minimum point in the kurtosis k{T^V), in 
Af{T, V) and in the third hexosis TLsiT, V): Kmin, ^ 3 ,mi™, 
Af*min 

(4) a singularity point in the pentosis n+{T, V),mU{T,V), 
in 0{T, V) and in the second hexosis 'H 2 {T, V): 7T+^oo, Woo, 
Ooo, '^ 2 , 00 - 

It is important to mention that the finite volume transi¬ 
tion point, using the connected Binder cumulant Bl{T, V), 
is given by the little maximum between the two 

minima. By against, the minimum of the Binder cumu¬ 
lant B 4 {T, V) : (i54)i„in obtained in [2111411140114211231124) . 
is just a particular point and not the true finite volume 
transition point. Obviously any particular point tends to 
the bulk transition point as V becomes large. The ap¬ 
parent discrepancy is completely due to the difference in 
the defintion of the Binder cumulant B 4 {T,V) and the 
connected Binder cumulant B%{T,V). The shift between 
{B 4 )^^^ and the true finite volume transition point in our 
model is different to those obtained by other models. This 
is probably due to the fact that our hpdf is very different 
from the double gaussian distribution used by Binder m 
and that considered in E3 ■ We therefore suspect that this 
shift is somewhere model dependent as quoted in El] ■ We 
will present a detailed study of this point in a forthcom¬ 
ing work. Finally, we can conclude that the finite volume 
transition point that appears as a particular point, the 
emergence of the linear correlation between different par¬ 
ticular points, and the possibility to use them to define a 
transition region, are the features of a universal behavior. 
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7 appendix 


From the general expression of the mean value (h") (T, V) 
(E3, we can easily deduce the first eight mean values : 

(h)(r, F) = (0)-/-oi (0) ^ (85) 


(h2)(T,F) = 


^01 ( 1 ) — -^01 ( 0 ) 

2 Lo3 (1) - 2 Lo3 (0) - 2Lo2 (0) - Loi(O) 


^01 (1) — ^01 (0) 


( 86 ) 


(h3)(r,F) = 
{h^)iT,V) = 
{h^){T,V) = 
(h6)(r,F) = 
{h^){T,V) = 


6 Lo4 (1) - 6 Lo4 (0) - 6 Lo3 (0) - 3Lo2 (0) - Loi(O) 


^01 (1) — ^01 (0) 


(87) 


24Lo5 (1) - 24Lo5 (0) - 24Lo 4 (0) - 12Lo3 (0) - 4Lo2 (0) - Loi(O) 

^01 (1) — ^01 (0) 

( 88 ) 

120Lo6 (1) - 120Lo6 (0) - 120Lo 5 (0) - 6OL04 (0) - 2OL03 (0) - 5Lo2 (0) - Loi(O) 

^01 (1) — ^01 (0) 

(89) 

72OL07 (1) - 72OL07 (0) - 72OL06 (0) - 360 £o 5 (0) - I2OL04 (0) - 30 Lo 3 (0) - 6£o2(0) - Loi(O) 

^01 (1) — Lqi (0) 

(90) 

7\rrr T/^ _ 5040Lo8 (1) - 5040Lo8 (0) - 5040Lo7 (0) - 2520Lo6 (0) - 840Lo5 (0) - 2 IOL 04 (0) - 42Lo3(0) - 7Lo2(0) - Loi(O) 


(91) 


^01 (1) — ^01 (0) 


(h«)(T,F) 


40320Lo9 (1) - 40320Lo9 (0) - 40320Lo8 (0) - 20160Lo7 (0) - 6720Lo6 (0) - I 68 OL 05 (0) - 336Lo4(0) 

^01 (1) — ^01 (0) 

56Lo3(0) + 8Lo2(0) + Loi(0) 

Lot (1) - Loi (0) ■ ^ 
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